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SUMMARY 


An  investigation  of  the  postbuckling  behavior  of  sandwich  plates  compressed 
beyond  the  general  instability  load  into  the  plastic  range  has  been  under¬ 
taken.  The  study  is  a  natural  extension  of  previous  work  which  led  to  the 
establishment  of  maximum  load-prediction  capability  relative  to  conventional 
plates.  The  purpose  of  the  present  investigation  is  to  assess  the  effects 
of  transverse  shear  deformations  on  the  maximum  strength  of  sandwich  plates 
when  the  primary  mode  of  initial  buckling  is  that  of  general  instability. 

The  analysis  utilizes  a  modified  version  of  Reissner's  variational  princi¬ 
ple  and  a  deformation  theory  of  plasticity  in  conjunction  with  computations 
carried  out  on  high-speed  computing  equipment.  The  results  indicate  that 
for  the  particular  plate  problem  investigated  (for  example,  aircraft  wing 
or  stabilizer  plate  element),  the  postbuckling  load-carrying  capability  is 
essentially  independent  of  transverse  shear  effects  in  the  presence  of 
inelastic  behavior  of  the  sandwich  faces.  The  implication  is  that  perhaps 
designers  will  not  have  to  be  as  conservative  in  sizing  sandwich  plate 
sections  for  maximum  load-carrying  capability  as  purely  elastic  theory 
results  would  indicate. 
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FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  a  continuing  effort  being 
undertaken  at  Stanford  University  for  the  U.S.  Army  Aviation  Materiel  Labo¬ 
ratories  under  Contract  DAAJ02-68-C-0035  (Task  LF162204A17002 )  to  establish 
accurate  theoretical  prediction  capability  for  the  static  and  dynamic 
behavior  of  aircraft  structural  components  utilizing  both  conventional  and 
unconventional  materials.  Predecessor  contracts  supported  investigations 
which  led,  in  part,  to  the  results  presented  in  References  2,3, 5,6,  and  8. 
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INTRODUCTION 


The  postbuckling  and  maximum  strength  analysis  of  square  and  rectangular 
plates  has  been  carried  out  by  Mayers  et  al.l>2,3  with  the  use  of  a  two- 
element  model  (see  Figure  l).  In  Reference  3,  Mayers  and  Nelson  state: 

"The  core  separating  the  faces  of  the  plate  is  considered  rigid  in  shear 
in  the  present  investigation,  but  provides  a  means  to  extend  the  current 
analysis  to  include  transverse  shear  deformation  when  postbuckling  behavior 
of  sandwich  plates  is  of  interest."  The  objective  and  contributions  of 
this  report  are  to  extend  the  earlier  work  and  to  establish  the  magnitude 
of  the  reduction  in  load-carrying  ability  of  postbuckled  sandwich  plates 
due  to  transverse  shear  effects.  The  study  of  sandwich  plate  and  shell 
behavior  has  been  undertaken  by  Mayers  et  al.  for  over  a  period  of  almost 
twenty  years  (References  4,5,  and  6).  With  ever-increasing  interest  being 
shown  by  designers  in  sandwich  construction  for  primary  structures  in 
aerospace  vehicles,  a  theoretical  prediction  capability  is  required. 

Although  much  progress  has  been  made  in  stability  analysis  of  sandwich 
structures,  very  little  information  is  available  for  the  prediction  of 
maximum  load-carrying  capability.  This  report,  then,  may  be  considered  as 
a  small  supplementary  portion  of  a  much  greater  research  effort  aimed  at 
developing  accurate  prediction  capability  for  the  bending,  buckling,  post¬ 
buckling,  and  maximum  strength  of  sandwich  plates  and  shells.  Most  of  the 
notations  used  here  are  the  same  as  those  in  Reference  3.  This  makes  it 
more  convenient  to  compare  the  equations,  expressions,  and  curves  presented 
here  with  those  appearing  in  References  1,  2,  3,  and  J.  These  comparisons 
show  the  effects  of  transverse  shear  on  the  maximum  load-carrying  capability 
of  infinitely  long  sandwich  plates  with  a  type  of  edge  condition  representa¬ 
tive  of  aircraft  wing  panels.  Although  calculations  have  been  carried  out 
only  for  very  long  panels,  the  theoretical  development  presented  herein  is 
established  for  retangular  plates  in  general.  Thus,  only  the  computational 
procedure  must  be  applied  to  produce  load-shortening  curves  for  such  plates. 
The  results  contained  herein  for  long  plates  indicate  that  transverse  shear 
effects  reduce  the  postbuckling  load-carrying  capability  to  much  the  same 
degree  of  magnitude  as  does  the  phenomenon  of  buckle  wavelength  change  in 
long  conventional  plates.  The  combination  of  postbuckling  waveform  modi¬ 
fication,  buckle  wavelength  change,  and  transverse  shear  flexibility  leaves 
little  postbuckling  load-carrying  capability  in  sandwich  plates  with  rela¬ 
tively  lightweight  cores.  For  aircraft  operating  in  the  low  and  medium 
subsonic  speed  range,  where  some  buckling  of  the  wing  surface  is  tolerable, 
only  a  more  detailed  comparison  study  can  disclose  which  of  the  two  types 
of  plates  (conventional  or  sandwich)  is  the  more  competitive  on  a  maximum- 
strength-to-weight  basis. 
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PROBLEM  DESCRIPTION  AND  GENERAL  THEORY 


STATEMENT  OF  THE  PROBLEM 


The  general  problem  considered  here  is  the  determination  of  the  buckling 
and  postbuckling  behavior  of  the  two-element  rectangular  sandwich  plate 
shown  in  Figure  1.  All  of  the  four  edges  are  simply  supported,  and  both 
the  loaded  and  unloaded  edges  are  constrained  to  remain  straight.  A  unit 
shortening  is  applied  in  the  x-direction,  which  up  to  buckling  corresponds 
to  a  uniform  axial  compression  Nx0.  For  convenience  and  with  justification 
obtained  from  the  quantitative  results  of  Reference  3,  the  translation  of 
each  point  in  the  y-direction  is  assumed  to  be  zero.  The  transverse  shear 
flexibility  of  the  core  in  the  xz-  and  yz-planes  is  considered,  thereby 
relaxing  the  rigid-core  assumption  inherent  in  References  1,  2,  and  3. 
Material  compressibility  is  neglected;  thus,  Poisson's  ratio  is  taken  as 
0.5  throughout  this  development. 

The  material  selected  is  2024-T3  aluminum  alloy  (see  Figure  2).  For  this 
aluminum  alloy,  the  material  constants  in  the  Ramberg-Osgood  uniaxial 
stress-strain  curve  representation 


«-f  ♦  «!>” 

are  K  =  3.14  x  1017 

n  =  8.60 

E  =  10.7  x  10^  psi 

with 

c  =  43.6  ksi. 

cy 


Two  different  plate  cross-section  geometries  are  studied.  For  the  purpose 
of  direct  comparison  with  earlier  work,  the  two-plate  cross  sections 
correspond  to  those  of  Reference  1,  which  are  designated  as 


and 


plate  1:  a  /a  =  0.6 
cr  cy 


plate  4 : 


a  /a 

cr  cy 


0.3. 


The  different  plate  critical  stresses  define  different  ratios  of  h/b.  For 
each  plate,  the  values  of  7  =  0,  0.25,  and  0.5  are  used  to  establish  the 
required  load-shortening  curves.  The  curves  for  7=0  (equivalent  to  a 
core  rigid  in  shear)  are  introduced  so  as  to  provide  a  check  with  the  load¬ 
shortening  curves  for  plates  1  and  4  (transverse  shear  effects  neglected) 
in  Reference  1.  The  main  jaurpose  of  this  research  is  to  study  the  variation 
of  (Nx)  /(Nx0)  with  u/ucr  in  the  postbuckling  range  to  assess  the 
magnitudeof  the  reduction  in  postbuckling  strength  caused  by  finite  shear 
flexibility  of  the  core. 


2 


I 


GENERAL  THEORY 

The  basic  problem  is  formulated  in  terms  of  a  variational  principle.  As 
demonstrated  in  Reference  7 ,  Reissner's  variational  principle  is  an 
alternative  to  the  potential  energy  principle.  In  References  3  and  8, 
Reissner's  principle  is  shown  to  be  a  better  alternative  than  the  potential 
energy  principle  for  nonlinear  problems  in  plates  and  shells.  In  terms  of 
both  stress  and  strain  components,  the  complete  Reissner  functional  for  the 
case  wherein  displacements  are  prescribed  at  the  boundaries  is  simply 


U"  =  /  {o  e  +a  e  +  t  y  +t  y  +t  y  -F'JdV 

J  x  x  y  y  xy  xy  xz'xz  yz 'yz 


where  F'  is  the  complementary  energy  density  for  a  nonlinear  elastic 
material  and  may  be  shown  to  be  given  by 
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€effd<Teff 


The  relation  between  effective  stress  a  __  and  effective  strain  e 

eff  eff 

relative  to  the  uniaxial  stress -strain  curve  for  the  material  may  be 
expressed  in  the  general  form  given  by  the  Ramberg -Osgood  formula 

,  _  !eff  ,feff  }n 

€eff  ~  E  + 

where  K  and  n  are  constants  which  depend  on  a  given  material.  With 
these  relations,  the  functional  can  further  be  expressed  as 


u"  =  /  (cr  e  +o  e  +t  y  +t  y  +t  y  -  [ 
J  x  x  y  y  xy'xy  xz'xz  yz'yz 


Teff  KE  ,aeffxn+111(,v 

+  n+rTT-'  ])dV  (2) 


Next,  the  coordinates  are  nondimensionalized  by  introducing  the  quantities 
£  =  2mx/L  and  r)  =  2ny/b;  the  total  volume  is  given  by  V  =  2t^Lb.  For 
convenience  in  computer  calculations,  the  functional  is  rewritten  as 

„„  r  r  <*'  <*'  T'  cr"  a"  t" 

U  =  /  /  {_£  6,  +  JL  e.  +  r  +  *  €"  +  JL  6"  +  JSL  y" 

EV  J  J  E  x  E  y  E  'xy  E  x  E  x  E  'xy 


n+1 


1  c  +  -v2  )  1r/  eff \  ,  f  eff\  i  K_r / 

£  Et  jcz  Vz  "  ^  E  't  E  V  2(n+l)  E  ' 


n+1 

+  (-|^)b  ]}didn  (3) 
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where  from  Appendix  I, 

■  i[(axl +  (ax)J 


.  A  1  2 
e'  =  u  +  -r  w 
x  ,x  2  ,x 


o'  =  k(a r  )  +  (a)  ] 
y  £  y  t  y  b 


.  .  l  2 

e'  =  v  +  t  w 
y  >y  2  ,y 


=  (t _ )  +  (O  ] 


xy  2  xy 


7'  =  u  +v  +ww 

xy  ,y  ,x  ,x  ,y 


»’x  -  •  («^l 

€x  ~  2^W,xx  '  \z,x^ 

e"  =  ~( w  -  y  •) 
y  2V  ,yy  yz,y 


7"  =  h[w  -  i(y  +7  )  ] 

xy  ,xy  2  xz,y  yz,x' 


and  the  functional  has  been  separated  into  the  contributions  of  the  direct 
and  bending  stresses  and  strains  in  the  two  faces,  and  the  shear  stresses 
and  strains  in  the  core. 

The  functional  given  by  equation  (3)  may  be  utilized  in  the  indirect 
variational  approach  to  develop  the  Euler  equations  and  boundary  conditions 
of  the  problem,  and  in  the  direct  variational  approach  (Rayleigh-Ritz 
procedure)  to  determine  the  generalized  stress  and  strain  distributions 
compatible  with  the  geometric  boundary  conditions  stipulated  for  the 
problem.  The  two  approaches  are  utilized  in  Appendixes  I  and  II,  respec¬ 
tively. 


k 


METHOD  OF  SOLUTION 


ELASTIC  SOLUTION 


In  the  elastic  case,  the  functional  (equation  (3))  may  be  written  as 
L  b 

U"  =  /  /  (N  (u  +  i  w^  )  +  N  (v  +  i  w2  )  +  N  (u  +  v  +  w  w  ) 
J  J  x'  ,x  2  ,x'  y'  ,y  2  >y/  xy'  ,y  ,x  ,x  ,y' 

0  0 

-M(w  -y  )-M(w  -y  ) 

x'  , xx  rxz,x'  y  ,yy  ryz,x' 

+  2M  [w  -  i(Y  +  Y  )]  -  rlr— [N^  +  N2  -  N  N  +  3N2  ] 
xy  ,xy  2  xz,y  'yz,x  4Et„  x  y 


x  y 


xy' 


. --^[M2  +  M2-MM  +3M2]+iGt(Y2  +  Y2  )}dxdy 

Et  h2  x  y  X  y  xy  2  c  cv  rxz  ryz/;  J 


00 


With  the  forcing  function  taken  as  the  unit  shortening,  u,  the  generalized 
displacements  and  forces  satisfying  at  least  geometric  boundary  conditions 
(see  Appendix  I)  are  assumed  to  be 

u  =  ux 

v  =  0 

w  =  C  sin  sin 

L  b 

Y  rnmx  .  iwy 
’xz  L  b 


Y  =  E  sin  — ^  cos  SSL 
yz  L  b 


(5) 


H  =  F  cos  +  11 

X  b  XO 


„  7;  2mnx 

N  =  G  cos  — — 

y  L 


N  =0 
xy 


y  j  nutx  .  rwy 
M  =  I  sin  -7—  sin  -r-*- 
x  Lb 
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M  =  J  sin  H22S  sin  HE 
y  Lb 

,,  77  mx  rmy 

M  =  K  cos  -7—  cos  — 

xy  Lb 


where  C,  D,  E,  F,  G,  I,  J  and  K  are  arbitrary  coefficients  whose  magnitudes 
are  to  be  determined  from  application  of  the  variational  principle.  Sub¬ 
stitution  of  the  assumed  functions  into  the  functional  (equation  (4))  and 
subsequent  integration  over  the  given  limits  yields  a  function  indepen¬ 
dent  of  x  and  y.  Application  of  Reissner's  principle  (variation  with 
respect  to  both  assumed  stress  and  strain  states)  yields  9  simultaneous 
equations  from  which  the  9  coefficients  can  be  determined  explicity  in 
terms  of  the  applied  unit  shortening.  A  detailed  development  is  given  in 
Appendix  II.  With  the  coefficients  established  for  any  given  value  of 
unit  shortening,  the  stresses  and  strains  at  each  point  of  the  plate  and, 
thus,  the  stress  distribution  are  known.  The  average  stress  or  load  per 
unit  length  at  any  cross  section  of  the  plate  may  be  plotted  versus  end 
shortening  to  produce  the  load -shortening  curve  in  the  elastic  postbuckling 
range.  Typical  load  distributions  (nondimensionalized  with  respect  to  the 
critical  load  (Nx0)  )  for  several  values  of  the  applied  shortening  param¬ 

eter  u  (nondimensionalized  with  respect  to  the  critical  shortening  ucr) 
are  shown_in  Figure  3.  The  average  load  ratio  (Nx)  /(Nx)  is  plotted 
against  u/ucr  in  Figure  4.  For  y  =  0  (infinite  sftear  rigidity  in 
core),  the  results  shown  in  Figure  4  are  identical  with  the  elastic - 
solution  results  obtained  in  References  1  through  3  for  square  or  infi¬ 
nitely  long  plates.  The  reduction  in  load-carrying  capability  for  y  = 

0.25>  0.5,  and  1.0  is  vividly  portrayed  in  Figure  4.  It  should  be  noted 
that  the  load-shortening  curves  appear  as  straight  lines  only  because  a 
one-term  approximation  has  been  utilized  to  represent  the  buckled  shape. 

In  reality,  this  one-term  approximation  is  exact  at  buckling  and  less  and 
less  accurate  as  u/ucr  increases.  As  shown  in  References  1  through  3, 
for  example,  the  straight-line  approximation  to  the  load -shortening  curve 
for  square  or  infinitely  long  plates  overestimates  the  exact  average  stress 
by  only  about  10  percent  when  u/ucr  is  as  great  as  9.  Since  the  large- 
deflection  plate  theory  is  being  utilized,  the  lateral  deflection  of  the 
plate  can  be  determined.  The  effect  of  transverse  shear  on  maximum 
deflection  is  shown  in  Figure  5. 

PLASTIC  SOLUTION 


In  the  plastic  range,  K  is  not  zero;  therefore,  the  exponential  terms 
will  appear  in  the  functional.  This  makes  it  impossible  to  apply  the 
variational  principle  to  obtain  an  analytic  solution  for  the  9  coefficients. 
However,  the  variation  can  be  carried  out  numerically.  For  a  certain  value 
of  unit  shortening,  approximate  numerical  values  of  the  coefficients  may  be 
obtained  by  taking  the  numerical  value  of  the  coefficients  of  the  elastic 
solution  as  initial  trial  values  for  the  Newton -Raphson  iteration  method  in 
the  plastic  solution.  When  the  differences  between  the  elastic  solution 
and  plastic  solution  values  of  the  corresponding  coefficients  become  large, 
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the  coefficient  values  of  the  preceding  plastic  solution  point  may  be  used 
as  initial  values  for  the  following  point.  The  principal  steps  are 
detailed  in  Appendix  III.  After  the  approximate  values  of  the  9  coeffi¬ 
cients  are  found  for  an  assigned  value  of  unit  shortening,  then,  using  the 
same  process  as  in  the  elastic  case,  the  stresses  and  strains  at  each 
point,  as  well  as  the  total  energy,  can  be  calculated.  Stress  distribu¬ 
tions,  average  stresses,  and  load -shortening  curves  follow  directly.  Load¬ 
shortening  curves  for  plates  1  and  k,  with  7=0,  0.25,  and  0.5  and 
plasticity  effects  included,  are  shown  in  Figures  6  and  7* 


RESULTS  AMD  DISCUSSION 


The  main  objective  of  this  report  is  to  establish  the  effect  of  transverse 
shear  flexibility  on  the  postbuckling  load-carrying  ability  of  both  elastic 
and  inelastic  plates.  As  can  be  seen  from  Figures  3  and  4,  based  upon  the 
all-elastic  analysis,  both  the  load  (stress)  distribution  across  the  plate 
and  the  average  load  are  markedly  changed  by  transverse  shear  effects.  In 
fact,  the  slope  of  the  load- shortening  curve  of  the  rigid  core  is  reduced 
from  l/2  (the  classical  result  for  conventional  plates  with  straight  un¬ 
loaded  edges)  by  more  than  20  percent  as  soon  as  7  reaches  the  value  of 
l/3.  Nevertheless,  Figure  4  indicates  that  a  substantial  load-carrying 
capacity  exists  in  the  elastic  sandwich  plate.  Now,  from  the  results  of  the 
inelastic  solution  (Figure  6  and  7)>  it  can  be  seen  that  plasticity  affects 
a  rigid-core  plate  to  a  much  greater  extent  than  it  does  a  plate  with 
finite  shear  flexibility  in  the  core.  In  other  words,  at  the  same  end¬ 
shortening  ratio,  the  effective  strain  levels  in  the  conventional  plate  are 
sufficiently  high  to  reduce  the  stresses  and,  hence,  the  average  stress, 
such  that  the  postbuckling -range  load  ratio  is  much  the  same  for  a  conven¬ 
tional  (7  =  0)  and  sandwich  plate  (7  =  0.25,  0.5).  The  results  shown 
in  Figures  6  and  7  indicate  that  the  usual  load-carrying  capability 
associated  with  conventional  plates  in  the  well-developed  inelastic  post- 
buckled  range  is  not  too  greatly  affected  by  the  presence  of  transverse 
shear  in  the  case  of  the  postbuckled  sandwich  plate. 

In  this  investigation,  considerations  of  the  face -wrinkling  (Reference  5) 
and  face-dimpling  (honeycomb  sandwich  only)  instability  modes  have  been 
neglected.  At  the  same  time,  both  waveform  changes  and  buckle  wavelength 
changes  during  postbuckling  have  been  neglected  (see  Reference  3).  The 
effects  considered  in  Reference  3  reduce  postbuckling  strength  and  cause 
the  occurrence  of  well-defined  maximum  loads.  Thus,  it  appears  from  the 
results  presented  herein  that  there  is  some  justification  for  designing 
sandwich  plate  structures  for  the  postbuckling  range,  since  the  plasticity 
and  waveform-wavelength  change  phenomena  (Reference  3)  reduce  load-carry¬ 
ing  capacity  (maximum  strength)  to  a  greater  extent  than  does  transverse 
shear.  Of  course,  in  any  design  exercise,  the  other  modes  of  initial 
buckling  must  be  considered. 


8 


CONCLUDING  REMARKS 


An  investigation  of  the  postbuckling  behavior  of  sandwich  plates  loaded 
into  the  plastic  range  has  been  undertaken.  In  the  absence  of  inelastic 
behavior,  it  is  demonstrated  that  transverse  shear  effects  significantly 
reduce  the  load-carrying  capability.  However,  when  inelastic  deformations 
are  taken  into  account,  the  load-carrying  capability  becomes  almost  inde¬ 
pendent  of  the  transverse  shear  flexibility  of  the  sandwich  core.  This 
result  indicates  that,  barring  initial  buckling  in  other  than  the  general 
instability  mode,  the  main  effect  of  transverse  shear  on  the  postbuckling 
load-carrying  capability  of  a  plate  is  to  precipitate  failure  (maximum 
strength)  at  a  lower  value  of  end  shortening,  with  little  change  in  the 
ratio  of  maximum  to  buckling  load.  Future  work  should  be  directed  toward 
establishing  the  magnitude  and  effect  of  inelastic  effects  in  the  core  on 
the  maximum  strength  of  plates. 
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eff  » 


Figure  2.  Compressive  Stress-Strain  Curve  and  Initial  Buckling 
Stresses  for  2024-T3  Aluminum-Alloy  Plates  Analyzed. 
(From  Reference  1,  Figure  2.) 


11 


?igure  5.  Variation  of  wmax/h  With  u/ucr  for  Several  Values  of  the 
Transverse  Shear  Stiffness  Parameter  y. 
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A,  A/;  y  =  0  (ref.  I  and  present  analysis) 

B,  B' :  y  =  0.25 

C,  C' :  y  =  0.50 


0  12  3  4 


_±L 

“cr 


Figure  6.  Variation  of  (Nx)bv/(Nxo)„„  With  u/ucr,  Plate  1>  cr  /o 

q  g  A  av'  cr  '  cr'  '  cr'  cy 
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APPENDIX  I 


EULER  EQUATIONS  AND  BOUNDARY  CONDITIONS  DERIVED  FROM 
THE  VARIATIONAL  PRINCIPLE 


1.  REISSNER  FUNCTIONAL 

For  a  sandwich  plate  as  shown  in  Figure  1,  the  Reissner  functionals  are 

U"  =  /  { (cr  e  )  +  (<j  €  )  +  (t  y  )  -  F !  ,  }dV 

t,b  J  X  x  t,b  y  y  t,b  *y  *y  t,b  t,b; 


for  the  face  plates  and 


U"  =  /  (t  y  +  x  y  -  F')dV 
c  J  xz  'xz  yz  [yz  c  ‘ 


for  the  core.  The  total  functional  is  given  by 

U"  =  U"  +  U"  +  U" 
t  b  c 

The  strain-displacement  relations  for  small  strain,  moderate-rotation 
plate  theory,  with  transverse  shear  effects  included  (see  Reference  7), 
become 


(e  ) 

»  *\ 

2 

w 

±  ^(w  -  y 

x  t,b 

,x  2 

2  ,xx  xz,x 

(e  ) 

V  +i 

2 

w 

±  ^(w  -  y 

y  t,b 

>y  2 

>y 

2V  ,yy  yz>y 

(y  ) 

=  U  +  V 

i-  w  w  ±  h[w 

*y  t,b 

,y  xy 

where  the  +  and  -  signs  correspond  to  the  top  and  bottom  faces,  respec¬ 
tively. 

The  complementary  energy  density  in  the  facings  and  core,  respectively,  is 
Pt,b  “  lE^x^b  +  (°A,b  ‘  (°x°y't,b  +  3(Txy\,b) 

k  •  sh-A +  A1 

c  J 
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With  these  relations,  the  functional  in  the  elastic  case  becomes 


U"  =  /{(<?)  [u  +  i  w2  +  — (w  -  y  )] 

J  x't  ,x  2  , x  2'  ,xx  rxz,x 

V 

+  (7Vty  *  Hy  +  t(u,yy  ’  W’ 

+  (o  )  [v  +  x  w2  -  x(w  -  y  )] 

y  b  >y  2  >y  2  >w  yz>y 

+  (t  )  [u  +  v  +ww  +  h(w  -  i(v  +  7  ))] 

xy  t  >Y  >x  ,x  ,y  ,xy  2wxz,y  yz,x"J 


+  (t  )  [u  +v  +  w  w  -  h(w  -  i(7  +  7  ))] 

xy  b  >y  >x  >x  >y  >xy  2 xz>y  yz>* 

-  +  3(Os 


y  t  x  t  y  t 


*y  t 


+  l*X  +  {\K  -  KW  + 

xb  y  b  x  b  y  b  xy  ^ 

+  t  7  +  T  7  -  ^t(t2  +  t2  )}dV 

xz  xz  yz  ryz  2G  xz  yz  ' 


Introduction  of  the  average  stresses  and  bending  moments 


o'  =  ±[(<J  )  +  (a)  ] 


x  2  x 


x  b 


a'  =  i[(cr  )  +  (a  )  ] 

y  2  y  t  y  b 

T1  =  ^[(t  )  +  (x  )  ] 

xy  2  '  xy't  xy'hJ 


M L  -  -tf  |t(dx)  -  (aj.  ] 


x'b' 


M 


-  -  KU 


y  t  y  b' 


M  =  t.  §[(t  )  -  (t  )  ] 

xy  f  2  xy  t  xy'bJ 


(6a) 
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into  equation  (6a)  yields 


U"  =  tfi/a<yx) t  +  (cJx}b](u,x  +  2  W,x)  +  [(ax)t  -  K\]  2(w,xx  ■  \z,x} 
00 

+  [(o  )  +  (a  )  ](v  +  ^  w2  )+[ (or  )  -  (o  )  ]-|(w  -  Y  ) 

y  t  y  b  >y  2  ,y'  y't  y  b  2V  ,yy  yz,y 

+  [(t  )  +  (t  )  ](u  +  V  +  w  w  ) 

xy  t  xy  b  ,y  ,x  ,x  ,y 

+  [(t  )  -  (t  )  ]h(w  -  ^[Y  +  Y  ]) 

'  xy'  '  xy  ,  v  ,xy  2L'xz,y  'yz,xJ/ 


4=[(o  )2  +  (o  f  -  (or  )  (o  )  +  3(x  )2 

2EL'  x't  y  t  x  t  y  t  xy't 


+  (of  +  (o  f  -  (a  )  (o  )  +  3(t  )2] 

b  y  b  byb  ^  b 

+  t— •[ T  Y  -  T  Y  -  -^r('f  +  f  )]dxdy 
tf  xzrxz  yz  yz  2Gc  xz  yz' 

or,  in  terms  of  local  average  stresses  and  moments, 


U"  =  t. 


(2o'(u  +  i  w2  )  +  2o'(v  +  i  w2  )  +  2t'  (u  +  v 

Xv  ,x  2  ,x  y'  ,y  2  ,y'  xyv  ,y  ,x 


i'i  i¥A 

+  w  w  )  -  t^(w  +  Y  )  -  r^(w  -  Y  ) 

>x  ,y'  tfv  ,xx  rxz,x'  tfv  ,yy  yz,y 

M 

+  2  ^[w  -  |(Y  +  Y  )]  -  i[d'2  +  O'2  -  O'O'  +  3t'2 

tf  ,xy  2  xz,y  'yz,x'  EL  x  y  xy  xy 

+  -  Jlf  +  M2  -  M  M  +  3M2  )] 

(th r  x  y  x y 


v  12  2 

+  t— [t  Y  +  t  Y  -  t^“(t  +  T  )])dxdy 

t  xz  rxz  yz'yz  2G  v  xz  yz,J1 

X 


2.  EULER  EQUATIONS  AMD  BOUNDARY  CONDITIONS 

The  vanishing  of  the  first  variation  of  equation  (6c)  with  respect  to  o', 

o',  t*  ,  M',  M',  M'  ,  u,  v,  w,  Y  and  Y  gives 
y*  xy  x  y  xy’  ’  ’  ’  ’xz  ’yz 
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t  T  t  T 

+  c  (y _ 2S5.)Bt  +  -£—(y _ ^)5t 

2t  ' 'xz  Go  '  xy  2t/  T®  G  Vz 

I  C  X  v- 

-(o'  +t'  )6u-(o'  +t'  )6v 

v  x,x  xy,y  y,y  xy,x' 

-[(o'w  )  +(o'w  )  +(t  w  )  +(r'  w  ) 

u  x  ,x',x  y  ,y  ,y  £y  ,y',x  xy  ,x  ,y 

+  +M  -2M  )]8w  +  rr-[-M  +M  +t  t 

04-  '  v  vv  nr  .inr  inr  .  vtr  '  J  9+.  L  t/-  .  V  YV  -  V  f»  ‘ 


2tfv  x,xx  y,yy  xy,xy 


2tf  y,y  xy,x  c  yz 


+  -f— [-M  +M  +t  t  ]8y  }dxdy 

2tf  x,x  xy,y  c  xz  xzJ 


d  l  ^  b 

2tf/{l,i6Ul  ^  +  2tf  /(VU|A)dX 


^tf/(oy6u|0,dx  +  2tf/(Vv|0,dy 

A  ^  A 


r  L 

•  /  {[2t„a'w  +2t„T'  w  +M  -M  ]B  |  }dy 

J  f  x  ,x  f  xy  ,y  x,x  xy,y  w  QJ 

0 

L  b 

■  /  ([2t„o'w  +2t„r'  w  +M  -M  ]5w  |  )dx 

J  u  f  y  »y  f  xy  ,x  y,y  xy,x 


L 

r 


■J  (Mx6YxzNdy  -  J 


u 

■/ 


+  /  {My&Yv7 |  )dx  -  r 

yz  0  J 


{M  6y  }dy 
xy  'yz  qj 


/Jj  r  u 

(M  6w  |  }dy  +  /  (M  8w  |  )dx 
*  j*  n  u  o 

A  V  A  W 
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(continued) 


» 


0 


(8) 


For  the  above  expression  to  vanish  for  arbitrary  states  of  stress  and 
displacement  consistent  with  the  geometric  boundary  conditions,  each  of 
the  terms  must  vanish  identically.  These  terms  lead  to  the  Euler 
equations  and  natural  boundary  conditions  as  follows : 


Stress -displacement  relations: 


A  1  2  dF'  A 

u  +  A  w  -  -r— p  =  0 

,x  2  ,x  do' 

’  x 


12 
v  +  7T  W 
>y  2  ,y 


3f' 


=  0 


u  +  V  +  w  w 
>y  ,x  ,y 


3F' 


=  0 


xy 


Moment -curvature  relations : 


2tf^W,xx"\z,x^  +  ° 

1  Sf1 

2t^w>yy  "  Yyz,y^  +  ST  "  ° 


i_[w  -  hy  +  y  )]  .  _ 0 

t  L  ,xy  2'7xz,y  7yz,x'J  3m 


(9) 


(10) 


Shear  stress -shear  strain  relations: 


t  T 


2tf'ryz  Gc  ' 


C-(Y 

Inplane  equilibrium 


x,x 


+  T1 


xy>y 


y>y 


+  T1 


xy,x 


-°  J 

equations : 


(11) 


(12) 
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Out-of -plane  equilibrium  equation: 

<y'w  +  2t'  w  +  cf'w  +  tt—(M 
x  > xx  xy  jxy  y  ,yy  2tf'  x,xx 

Core  shear  equilibrium  equations: 


-M  +  M  +  t  t  =0 

x,x  xy,y  c  xz 

-M  +  M  +  t  T  =0 

y>y  xy>x  c  yz 


Stress  (natural)  boundary  conditions: 


2M  +  M  )  =  0 
xy,xy  y,yy 


t' 

xy 

=  0 

at 

<< 

ii 

o 

\a 

y 

=  b 

M 

X 

=  0 

at 

O 

II 

K 

X 

=  L 

M 

y 

=  0 

at 

II 

o 

y 

=  b 

The  geometric  boundary  conditions  that  have  been  enforced  are 


w  =  0 
w  =  0 
u  =  0 
u  =  uL 
v  =  0 
v  =  0 


x  =  0,  x  =  L 

y  =  0,  y  =  b 

x  =  0 

x  =  L 

x  =  0,  x  =  L 

y  =  0,  y  =  b 
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The  variation  of  U"  with  respect  to  the  parameters  gives 
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#_  u"  =  Lbu  +  £  ^  ^  -  ssr  H 

N  o  L  2Et_  xo 

xo  f 


(l8a) 

(18b) 

(l8c) 

(l8d) 

(l8e) 

(l8f) 

(l8g) 

(l8h) 

(l8i) 


First,  consideration  is  given  to  the  linear  (excluding  N  )  part  of 
equation  (l8a)  and  equations  (l8b),  (l8c),  (l8f),  (l8g)  and  (l8h).  These 
6  equations  are  linear  and  homogeneous  with  respect  to  the  parameters 
C,  D,  E,  I,  J  and  K.  A  nontrivial  solution  exists  only  if 
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2.  CRITICAL  STRESS 

With  the  introduction  of  the  parameters 
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expansion  of  the  above  determinant  leads  to 
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or 


“x  -  -  -a  <2°> 

x  3[l+r27(r+l)] 

For  the  critical  value,  n  =  1  and  k  is  minimized  with  respect  to  3,' 
that  is,  dk^/dP  =  0.  This  gives 


Substitution  of  the  initial  wavelength  parameter  value  into  equation  (20 ) 
gives 


(k  ) 

'  x  cr 


_ 4 _ 


(21a) 


If  y  =  0,  that  is,  approaches  infinity,  then 

(k  )  I  =  -4  (21b) 

x  cr  ly-o 

The  negative  sign  simply  indicates  that  u  has  been  taken  positive  in 
tension. 


3.  STRESS  DISTRIBUTION 

The  values  of  all  arbitrary  parameters  may  be  found  in  terms  of  C  as 
follows:  First,  from  equations  (l8b),  (l8c),  (l8f),  (l8g)  and  (l8h), 
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Second,  equations  (l8d),  (l8e)  and  (l8i)give 
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Substitution  of  P,  G,  I,  J  and  K  into  equation  (l8a)  then  yields 
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anc^  n  =  1,  the  critical  value  for  k  exists  and 
1 1_Y  x 


(?)  -^Quk)  -k) 


3(1+7  )  x  cr 


32(1-7) 


3(l+72)(l-*7) 


KH 

cr 


(25) 


Here,  it  should  be  mentioned  thax  when  equation  (25)  is  substituted  into 
the  expressions  for  the  parameter  coefficients,  all  P's  should  take  the 
value  of  v(l+7)/(l-7)  when  the  plate  is  long.  For  finite  length  plates, 
integer  values  of  p  must  be  utilized,  just  as  in  the  case  of  conventional 
plates,  to  find  the  lowest  critical  k^.  Thus, 

N  =  F  cos(%£)  +  N 
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N  /.  ,  \2  r  N  1  _  N 
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To  find  the  relation  between  u  and  N**  N  must  be  substituted  into 
equation  (l8a)  also.  Then,  equation  (18a)  becomes 
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This  leads  to  the  critical  value  of  u, 
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when  C  is  set  to  zero,  n  =  1,  and_  3  =  s/(l+y)/(l-y).  It  can  be  shown 

here  that  this  is  in  agreement  with  u  =  cr  /E  where  a  =  (N  )  /2t„. 

°  cr  cr'  cr  xo'cr'  f 

Now,  C,  in  terms  of  u  and  u  ,  is 
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On  the  other  hand,  equation  (l8i)  may  be  rewritten  as 
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Next,  N  can  be  expressed  in  terms  of  u  and  u  as 
x  cr 

N  =  F  cos  (“2Z)  +  N 
x  '  b  xo 
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For  y  =  0,0.25,0. 5,1.0,  the  relation  between  (N  )  /(N  )  and  u/u 

is  shown  graphically  in  Figure  4.  *  i: 

4.  MAXIMUM  DEFLECTION 

From  equation  (5);  the  lateral  deflection  is 

w  =  C  sin  ^  sin  SSL 
L  D 


The  maximum  value  occurs  at  x  =  L/2m  and  y  =  b/2.  Thus, 


From  equation  (34),  w  /h  can  be  expressed  in  terms  of  u/u  ,  that 
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Equation  (38)  agrees  with  equation  (25)  when  C/h  is  expressed  in  terms 

of  u/u  . 

'  cr 

For  different  values  of  7(0, 0.25,0. 5,1.0),  the  variation  of  w  /h 
with  u/ucr  is  shown  in  Figure  5-  max 


APPENDIX  III 


PLASTIC  ANALYSIS  AND  NEWTON -RAPHSON  ITERATIVE  TECHNIQUE 

In  the  plastic  case,  K  is  not  zero.  Therefore,  analytic  expressions 
for  the  assumed  displacement  and  stress  coefficients  are  not  obtainable. 
The  Newton-Raphson  iterative  technique  is  utilized  to  obtain  the  numerical 
solutions  corresponding  to  vanishing  of  the  first  variation  of  the  func¬ 
tional  given  by  equation  (3). 

In  the  computer  program,  it  is  more  convenient  to  write  the  functional 
corresponding  to  equation  (6)  for  the  elastic  case  in  the  following  form: 


Programming  of  this  equation  constitutes  the  main  subroutine.  Also,  the 
following  expressions  are  used  in  the  program.  New  coefficients, 
designated  by  a  tilde  are  introduced.  For  the  elastic  portion, 
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It  should  be  noticed  that  p's 

here 

must 

take  the 

va_ue  of  nT(  1+y)/(1-'y) 

so  that  the  coefficients  will  agree  with  those  taken  earlier  in  the 
expressions  (25)  and  (29 ) •  This  is  also  true  for  the  following  expressions. 
Then,  the  stresses  and  strains  in  equation  (39)  are 
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Because  of  symmetry,  the  functional  integration  is  performed  for  a 
quarter  of  the  plate.  For  a  specified  material,  for  example,  one  of 
2024 -T3  aluminum  alloy  (Figure  5), 

K  =  3.14  x  1017 
n  =8.6 

E  =  10.7  x  106  psi 

a  =  43.6  ksi 

cy 

For  plate  No.  1, 
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whereas  for  plate  No.  4, 
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From  these  given  values,  uc^  and  h/b  are 
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h  /  5(,cr  Hy  /5<T cr 

5  '  /A(k  )  =  2it  J  E 

1  '  x  cr 

Hence,  with  assigned  values  of  7  and  u/u  (elastic  case),  the 
stresses  and  strains  at  each  point  can  be  determined,  and  the  integration 
of  equation  (39)  can  be  made  numerically  by  applying  Simpson's  rule  in 
two  dimensions.  The  corresponding  coefficients  of  the  plastic  solution 
can  be  found  numerically  by  using  the  elastic-solution  coefficients  (or 
known  plastic  coefficients  of  a  previously  calculated  neighboring  point) 
as  initial  values  and  by  applying  the  Newton-Raphson  iterative  technique 
to  obtain  a  stationary  value  of  the  functional.  This  involves  the 
solution  of  9  simultaneous  nonhomogene ous  linear  equations. 

The  functional  U"/EV  after  integration  is  de^igned^by^  U",  ^and  the 
9  coefficients  of  the  plastic  solution  ((5^  15,  E,  ?,  "5,  J,  K,  and  iJ  ) 

are  designated  by  Xj_ .  The  quantities  dU"/dx^  and  c^2U"/dxi^x.  are° 
designated  by  and  Uj'^j .  Then,  the  set  of  equations  is 

u;,11ax1  +  u;'12ax2  + . +  u;19ax9  =  -  u;^ 

u;21ax2  +  u;22ax2  + . +  u;29ax9  =  -  u;'2 

u;9iAx1  +  u;92ax2  + . +  uv99ax9  =  -  u;9 

where  Ax^  are  the  unknown  increments  in  the  iterative  process,  and  all 
of  the  U"  ,  U"  can  be  found  by  the  numerical  expressions 
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"  (8x.)2 

U"(x.  +  5x.,  x.  +  8x.)  -  U"(x.  +  8x . )  -  U"(x,  +  8x.)  +  U" 

Tin  _ 1 _ 1 _ £ _ J _ 1 _ _ jJ _ J _ 

,ij  6xi5xj 

-3  —  —  — 

The  increments  5x^  are  taken  as  x.^  «  10  .  All  of  the  U",  Uj_',  are 

evaluated  with  the  assigned  initial  coefficients.  This  involves  signif¬ 
icant  computational  effort  since  different  values  of  7  for  the  desired 
range  of  u/ucr  have  to  be  assigned. 
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After  sufficient  iterations,  the  approximate  numerical  value  of  each 
coefficient  can  be  found  to  a  desired  degree  of  convergence.  Then,  the 
stresses  and  strains  at  each  point  of  the  plate  can  be  calculated.  Thus, 
the  average  value  of  <jx  can  be  determined  by  numerical  integration  of 
the  stress  distribution  across  y.  The  variation  of  (ax)aY/(ax) 

(which  in  this  problem  is_the  same  as  (Nx)  /(Nxq)  )  with  respect  to 
the  assigned  value  of  u/ucr  can  be  plottetvin  graphical  form.  This 
curve  determines  the  load-carrying  capability  of  the  sandwich  plate  in 
the  postbuckling  range.  Load -shortening  curves  for  plates  1  and  4,  with 
plasticity  effects  included,  are  given  in  Figures  6  and  J,  respectively. 

With  reference  to  the  actual  computer  programming,  it  was  found  that  the 
more  digits  employed  (in  Algol  W,  long  real  numbers  must  be  used),  the 
better  the  convergence.  The  Newton-Raphson  iteration  technique  is 
basically  a  very  effective  method  of  numerical  minimization.  However, 
initially,  seven-digit  calculations  led  to  a  very  poor  convergence.  A 
simple  test  program  based  on  a  quadratic  form  in  three  variables  showed 
that  the  convergence  became  remarkably  rapid  when  15  digits  were  employed. 
Of  course,  since  the  actual  problem  reflects  9  arbitrary  coefficients  and 
a  high  degree  of  nonlinearity  in  the  functional  that  was  varied,  con¬ 
vergence  was  slower  but,  nevertheless,  quite  satisfactory. 
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APPENDIX  IV 


EFFECTIVE  STRESS  AND  STRAIN 

The  complementary  energy  density  for  linear  elastic  behavior  of  plates 
is  (see,  for  example,  Reference  3) 

F'  =  4=r{(X2  +  cr2  -  2pcr  a  +  2(1  +  m->Y2  ) 

2E  x  y  x  y  v  rxyJ 

This  can  be  easily  verified  by  enforcing  the  generalized  stress-strain 
relations 
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(42) 


Then,  from  any  one  of  References  1  through  3, 
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(43) 


In  linear  elasticity,  aeff/eeff  =  E>  thus, 
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Integration  gives 


tf2  =  F'  =  Jtt{<j2  +  a2  -  2p0  cr  +  2(1  +  p)y2  } 
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(45) 


In  nonlinear  elasticity  (References  1  through  3),  the  same  equations 
hold.  For  p  =  0.5,  the  expression  for  the  effective  stress  is 


eff 


2  2 

=  .  a  +  a  -  a  a  +  3y 
x  y  x  y  fxy 


(46) 


This  expression  may  be  transformed  into  effective  strain  simply  by  using 

the  relations  (42),  with  E  replaced  by  E  ;  that  is, 

s 
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